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Abstract. We apply the stochastic approach to the calculation of the Reggeon Field Theory (RFT) 
elastic amplitude and its single diffractive cut. The results for the total, elastic and single difractive 
cross sections with account of all Pomeron loops are obtained. 
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THE STOCHASTIC APPROACH 

The elastic scattering amplitude of the RFT [1] is given in terms of quasiparticle ex- 
changes (Pomerons and Reggeons), its cuts determine cross sections of various inelastic 
processes. While the lower energy data on total and elastic cross sections can be de- 
scribed in terms of non-interacting Pomerons and Reggeons, the high energy behaviour 
of total, elastic and high-mass diffractive cross sections makes the account of Pomeron 
interactions and loops in the RFT fits absolutely essential. For this task we apply a 
method which is based on the approach known as reaction-diffusion (RD) or stochastic. 

It was observed [2] that a stochastic system of classical particles in the 2-dimensional 
plane admits a field-theoretical description with the Lagrangian of the Gribov-Regge 
theory with triple and 2^2 interaction terms: 

^ = - - a^(Vb0t)(Vb0) + A(/)t(/) + /r3p(/)^0(0^ + (/>) +Z0^'</>' (1) 

The classical partons are allowed to move chaotically (characterized by diffusion coef- 
ficient D), split, A ^ A + A, with a splitting probability per unit time A, or die, A ^ 0, 
with a death probability mi. When two partons are brought within the reaction range a 
due to the diffusion, they can pairwise fuse, A + A ^ A, or annihilate, A + A ^ with 
the rates v and m2 correspondingly. A set of inclusive ^'-parton distributions: 

f^{y;^^) = Y,n^f^Jd^NPN{y',^^^ (2) 
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obeys the same evolution equations as the set of the exact Green functions for the La- 
grangian (1). Here p^(y\^^) are the symmetrized probabihty densities with normahza- 
tion Y^NwJ d^NpN{y\^N) = lLMPN{y) = l and ^ = {zi, . . . ,z^}. Evolution time for 
partonic system in RD approach is dual to the rapidity in the RFT. 

Phenomenological parameters of the Lagrangian (1) have direct correspondence with 
the rates of the stochastic system (see table 1). The parton interaction distance a serves 
as a regularization parameter for the Pomeron loops. For given values of the coupling 
r^p and the scale £ = na^ the quartic coupling x can be varied. 



TABLE 1. Relation between the parameters of the RFT and those of the stochastic approach. 



RFT a' 


A 


r^p, P splitting vertex 


r^p, P fusion vertex 


X, 2P 2P 


RD-approach 


D 


A — mi 




(m2 + ^v)v/e 


^(m2 + v)£ 



This equivalence allows to compute the exact Green functions of the RFT with the 
account of all loop contributions following a Monte-Carlo evolution of the RD system. 



AMPLITUDE CALCULATION 

The paper [3] suggested a straightforward approach to the amplitude calculation. One 
has to compute the ^-point Green functions in convolution with hadron-fz-Pomeron 
vertices (grey blocks in fig. la) for the projectile (fs) and target(/^) and perform their 
numerical convolution according to the general rules of the Reggeon field theory at some 
linkage point y in rapidity. In particular, the elastic scattering amplitude is : 

Here g are some narrow functions normalized to / g{b)(fih — £. 

The most efficient way to compute the amplitude is its computation on the event by 
event basis with subsequent Monte-Carlo average setting the linkage point to the target 
rapidity. The inclusive fz-parton distributions at projectile and target rapidities coincide 
with the hadron-/z-Pomeron vertices [4]: 

fs{y = 0;^s)^^isPs{^s) = £'^^^^'\^s). (4) 

This gives the distribution of partons at zero evolution time in number and positions in 
the transverse plane. In accordance with it the initial configurations of partons should be 
generated within the first step of the Monte-Carlo averaging procedure. 

Upon the Monte-Carlo evolution of initial random parton configuration one gets a 
set of N partons at certain positions bi in the transverse plane. The event realization of 
inclusive distribution is thus fs^^^\^s) = L{/i.../Jg{i...a^} ^(^i ~ ^/i) • • • ^{^s — ^/J and 
upon convolution with the set of target-^-Pomeron vertexes leads to: 

7;lpie(b) = i:(-ir^A.e^ I Ps{\-h,...,±i-h). (5) 

s=\ ii<i2...<is 

The actual value of the elastic amplitude as a function of the impact parameter b is 
computed by making Monte-Carlo average of (5). 




FIGURE 1. Diagrams taken into account in calculation of elastic amplitude (a-c) and its SD cut (d-g). 

A single diffractive cut of the amplitude for a given value of the rapidity gap _ygap (with 
elastic contribution included) can be computed as a sum of two terms ^: 

rSD-t(b,F,3;g,p) = 2r\hj) - r'(b,F,ygap). (6) 

The term T\h, Jgap) is computed in the same way as the elastic amplitude by making 
a Monte-Carlo average of (5) with only the distinction in preparation of the projectile- 
associated set of partons. The evolution starts with two sets which evolve independently 
up to the evolution time jgap corresponding to the width of the rapidity gap. At that point 
the resulting partons are combined into a single set which further evolves in the standard 
way from y^^^ up to the target rapidity Y . 



PARAMETERS AND RESULTS 



We use the approach described above to compute the elastic amplitude in impact param- 
eter representation and its single diffractive cut with the account of Pomeron trajectories 
only. For the parameterization of the proton-fz-Pomeron vertexes we use the two-channel 
eikonal approximation with the same values of the channel weights C\ —C2 =0.5 and 
of the parameter r; = 0.55 for the relation between the channel couplings to the Pomeron 
j8i/2 = ( 1 =t ^ ) j8o as in [5] . We however use a Gaussian parametrization of these vertexes 

and neglect the real part of the Pomeron exchange amplitude, \mAp{b) — T^^{b). The 
triple coupling value = 0.087 GeV~^ is taken according to [6]^. 

In order to obtain a better description of 
the data at lower energies we add a lowest 
order contribution of secondary trajecto- 
ries of positive and negative signature (see 
fig. 1). This contribution is added numer- 
ically to the all-loop Pomeron exchange 
amplitude given by the Monte-Carlo com- 
putation as described above making use of 
the Regge factorization. We also keep the 



Trajectory 


P /?+ R- 


Ao/+/- 
a', GeV-2 
R^, GeV-2 
/3o/+/-,GeV-i 


0.195 -0.34 -0.55 
0.154 0.70 1.0 
3.62 3.0 5.2 
4.7 4.05 2.59 



TABLE 2. Parameters of the trajectories and their 
couplings to proton. 



real part for the secondary Reggeon contribution. 

We perform a fit of total and elastic cross sections together with elastic scatter- 
ing slope. For this we in addition fix apriori the regularization scale a = 0.036 fm= 
0.182 GeV~^ and 2 ^ 2 coupling x by setting v = 1 /2A and m2 = (see tab. 1). This 
gives X — 2.87 x 10~^ GeV~^. This set is referred to as "set 3" in fig. 2. The parameters 
of Regge trajectories used and couplings to proton are listed in tab. 2. 

To illustrate the dependence of the result on the regularization scale and the quartic 
coupling we perform calculation also for the other parameter sets which differ from 



We are currently preparing a paper with a complete derivation of eq. (6). 
Our normalization of (see [4]) differs from that in [5] and [6]. 
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FIGURE 2. Total cross section, elastic cross section and elastic scattering slope. 
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FIGURE 3. Left: single diffractive cross section for < 0.05^ (jgap = 3, upper curves) and low-mass 
(ygap = 2ln^/s, lower curves) dissociation. Center, right: unintegrated profiles for inelastic and single 
diffractive cross sections, high- and low-mass (for set 3). 



the original (set 3) by values of a and x ^i^d have the same parameters of the Regge 
trajectories (sets 1, 2, 4). Namely we take (^1=^2 = 0.091 GeV~\ = 0.182 GeV~^; 
Xi=Xa = 1.435 X 10-2 GeV-2, X2 = 1.11 x 10"^ GeV-^. 

In fig. 3 (left plot) we show results of the diffractive cross sections calculation for the 
sets 1-4 (for the accounted graphs see fig. 1 d-g). Plots for the inelastic and diffractive 
profiles indicate that at high energies the inelastic diffractive contribution comes domi- 
nantly from the ring on the periphery of the disc. This is consistent with the logarithmic 
growth of the all-loop inelastic diffractive cross sections with collision energy. 
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